 use "$data\women_analysis.dta", clear

 global controlwoman mbase_ageyears  fdstotal rantime ranerror numberkids
 global controlhh land members adults mbase_electrified  mbase_OB mbase_scheduled
 
 estimates clear
  foreach i in  Pleave2 Pmarketin2 Pmarketout2 Phealth2 Pnatal2 Pcallnatal2 Pgroup2 Pmela2{
  ivregress 2sls `i'  (participatedTA_IV=ITTanalysis) i.strata $controlwoman $controlhh if sequence==1, vce(robust)
  estimates store `i'
  su `i' if sequence==1 & ITTanalysis==0
  estadd scalar controlmean=`r(mean)'
 }
 
 esttab * using "$results/TableB10.xls", stat(controlmean N r2) keep (participatedTA_IV) b(3) se(3) replace
 
 ***FDR q-values
 *The FDR Q-values are computed following  Benjamini and Hochberg (1995) within each category.
 
 exit
 
 





